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Procede de determination des contraintes, 
deformations, endommagement de pieces constitutes d'un 

materiau solide 



La presente invention concerne un procede pour 
determiner les effets d'une sollicitation mecanique sur 

5 un objet constitue d ! un materiau solide. 

Dans de nombreux domaines techniques , on a besoin de 
determiner l f incidence de sollicitations mecaniques sur 
des pieces , notamment des pieces metalliques, soit afin 
d'optimiser des procedes de fabrication comportant des 

10 operations de mise en forme, soit afin de prevoir le 
comportement en service de pieces dans des conditions 
complexes . 

Les procedes de mise en forme concernes sont les 
procedes de mise en forme par deformation plastique tels 
15 que l'emboutissage de tales, I'etirage de tubes, le 
forgeage de pieces de mecanique. Ces procedes peuvent 
avoir une incidence importante sur les proprietes du 
materiau constitutif des pieces concerntes, et il est 
important de pouvoir prevoir cette incidence de fagon a 
20 permettre d'evaluer a priori le comportement en service 
des pieces obtenues. 

Les sollicitations en service des pieces sont 
diverses et correspondent & toutes les sollicitations que 
l f on peut rencontrer, par exemple a l'effet d 1 efforts 
25 permanents, l'effet d 1 efforts appliques progress ivement a 
des vitesses plus ou moins elevees, d f efforts appliques 
de fa<pon periodique, d f efforts appliques dans des 
conditions de temperature ambiante ou bien d 1 efforts 
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appliques dans des conditions de temperature soit 61evee 
soit tres basse. 

Afin d'evaluer les effets de ces sollicitations 
mecaniques sur des pieces, on utilise de fa<?on connue des 
5 methodes de calcul num^rique, en mode elastique ou 
an61astique. Ces methodes de calcul sont, notamment, les 
methodes par differences finies et les methodes de type 
Galerkin reposant sur une decomposition sur des fonctions 
de base dont un type particulier est la methode par 
10 elements finis. Par "mode anelastique" , on entend tout 
mode de comportement qui n'est pas elastique reversible, 
c'est a dire tout mode de comportement plastique, 
plastique dependant du temps, ou viscoplastique, 
plastique dependant de la temperature, plastique 
15 resultant de transformation de phases. Ces calculs 
utilisent une modelisation du comportement du materiau 
sous l'effet des contraintes, representee par une loi 
reliant les contraintes aux deformations. De fagon 
classique, on utilise des lois de comportement 
20 macroscopiques qui correspondent a des representations 
phenomenologiques des phenomenes de deformation, c f est-a- 
dire des representations qui ne tiennent pas compte des 
phenomenes physiques fondamentaux sous- jacents . De tels 
modeles ont 1 1 avantage d'utiliser un nombre de variables 
25 relativement limite ce qui conduit a des calculs 
suffisamment peu volumineux pour pouvoir etre utilises 
dans les processus de conception industriels. En 
revanche, ce type de modele presente 1 1 inconvenient de ne 
pas permettre de modeliser convenablement l'effet de la 
30 microstructure des materiaux sur leurs propriety, et en 
particulier sur les proprietes anisotropes, et plus 



WO 2005/015167 



3 



PCT/FR2004/002012 



encore sur 1 1 Evolution de ces propriety anisotropes sous 
.l ! effet des sollicitations et des deformations. 

Afin de remedier k cet inconvenient , il a et6 propose 
d ! utiliser des modeles de comportement de materiaux, dits 
5 "modules polycristallins" , dans lesquels le materiau est 
represents par une pluralite de grains ayant la structure 
cristallographique et les propri6t6s d'un grain du 
materiau. Chaque grain est caracterisS notamment par son 
orientation et par les systemes de glissement propres a 
10 sa structure cristallographique. Dans ces conditions, le 
nombre de systemes de glissements est important puisqu'il 
va de douze pour un cristal cubique k faces centrees a 
trente pour un cristal hexagonal compact. De tels modeles 
ont l'avantage de permettre une representation fine du 
15 comportement du materiau, mais 1 1 inconvenient de 
nScessiter 1 1 utilisation d f un nombre tres important de 
variables. Ce grand nombre de variables conduit k des 
temps de calculs absolument prohibitifs qui ne permettent 
pas 1 1 utilisation de tels modeles dans un contexte 
20 industriel. Enfin aucun de ces modeles ne permet de 
prendre en compte les endommagements microscopiques des 
materiaux, et done ces modeles ne permettent pas de 
prevoir les ruptures en service. 

Le but de la presente invention est de rem6dier a ces 
25 inconvenients en proposant un procede pour determiner les 
effets d'une sollicitation mecanique sur un objet 
constitue d'un materiau solide qui permette de simuler de 
fagon fine devolution des proprietes d'un materiau 
anisotrope, en particulier d'un materiau polycristallin, 
30 en utilisant un nombre suffisamment restreint de 
variables pour que les temps de calcul soient 
raisonnables et permettent de faire des simulations dans 
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un contexte industriel, tout en permettant de prendre en 
compte les endommagements du materiau. 

A cet effet, 1' invention a pour objet un procede pour 
determiner les effets des sollicitations mecaniques sur 
5 une piece constitute d f un materiau solide, selon lequel 
on calcule les deformations et les contraintes engendrees 
par la sollicitation mecanique en une pluralite de points 
de la piece au moyen d'une methode de calcul numerique en 
mode anelastique, dans lequel le comportement du materiau 
10 solide est reprtsente par un module de comportement 
microscopique polycristallin utilisant une pluralite de 
briques "grain" ayant chacune une orientation et une 
fraction volumique, dont les deformations sont 
determinees a partir d'une pluralite de systemes de 
15 glissement propres au materiau solide, la trace du 
tenseur des deformations plastiques microscopiques de 
chaque brique "grain" etant nulle, les orientations 
distinctes et les fractions volumiques correspondant a 
des briques "grain" etant representatives de 
20 1 1 anisotropie du materiau solide. Selon l 1 invention, le 
nombre de briques "grains" est inferieur ou 6gal a dix et 
le nombre de systemes de glissement utilises dans le 
modele est inferieur ou egal a six. En outre, le modele 
microscopique polycristallin peut eventuellement utiliser 
25 jusqu'a sept briques "joints de grains" ayant chacune une 
orientation et une fraction volumique, et dont les 
deformations sont determinees a partir d'au moins un 
systeme d'ouverture, la trace du tenseur des deformations 
plastiques microscopiques de chaque brique "joint de 
30 grain" etant positive. Le modele peut egalement utiliser 
eventuellement une brique "trou" ayant une fraction 
volumique variable, dont les deformations plastiques 
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sont purement volumiques et dont la trace du tenseur des 
deformations plastiques microscopiques est positive. 

De preference, les six systemes de glissement 
associes k une brique "grain" definis dans un rep£re 
5 orthonorme associe a la brique "grain" definis par les 
normales au plan de glissement et par les directions de 
glissement sont tels que le double de la matrice 
d f orientation associee & chaque systeme de glissement 
n'est constitue que de "0", de "-1" et de "1". 
10 Pour un materiau isotrope a structure cubique les 

deformations sont determinees de preference a partir de 
deux families de trois systemes de glissement, une 
premiere famille correspondant aux faces du cube 
cristallin, et la deuxi^me famille correspondant aux 
15 plans de 45° de ces faces. De preference le nombre de 
briques "grain" est 6gal a sept et les briques "grain" 
sont obtenues par symetrisation de deux briques "grain" 
primitives situees sur les grands cercles a 45° de la 
figure de poles. 

20 On peut en outre determiner devolution de la texture 

du materiau en calculant la rotation Q du r6seau 

cristallin de chaque brique "grain" par rapport a un 
repere corotationnel. Lorsque le materiau est isotrope a 
structure cubique, le nombre de briques "grain" est alors 

25 de preference 6gal a dix. 

On peut egalement determiner 1 ' avancement d'une 
transformation de phase. Pour cela, on peut attribuer a 
chaque brique "grain", des caracteristiques fonction de 
l 1 avancement de la transformation de phase. 

30 D'une fagon generale, les orientations et les 

fractions volumiques des briques "grain" peuvent etre 
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determinees par ajustement a partir d'essais de traction 
biaxiale et de traction torsion. 

En l 1 absence de brique "trou", le modele de 
localisation pour le passage du niveau macroscopique au 
5 niveau microscopique pour une brique "grain" , permettant 
de calculer le tenseur des contraintes microscopiques a 

— 8 

dans la brique "grain" en fonction du tenseur des 
contraintes macros copiques E et du tenseur £ g p des 
deformations plastiques microscopiques dans la brique 
10 "grain" et du tenseur B P des deformations plastiques 
macroscopiques s'6crit : 

g =E+2n(l-p) [rid® l)+otl] (E p .e p ) 

r=g =3 1 S3 IS = — — g 

Dans cette formule cc et ti sont deux coefficients 
15 dependant d'une part du materiau et d 1 autre part de la 
valeur maximale atteinte au cours de la deformation par 
le deuxieme invariant de von Mises de la deformation 
plastique macroscopique, 

De preference, lorsque les calculs sont des calculs 
20 par 616ments finis, on associe une seule brique "grain" 
et eventuellement une brique "trou" ou une brique "joint 
de grain" a chacun des points d' integration de la m^thode 
de calcul par Elements finis ; la repartition des briques 
entre les points d 1 integration de la methode de calcul 
25 etant effectu^e de telle sorte que cette repartition soit 
homogene . 

La methode de calcul par Elements finis utilisee peut 
etre une methode de calcul mecanique ou thermo-mecanique 
faisant intervenir le temps. 
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Dans ce cas, lorsqu'on utilise une brique "trou", la 
fraction volumique cie la brique "trou" n'est pas nulle a 
1 ! instant initial des calculs. 

Lorsque le nombre de briques "joint de grain" ou de 
5 briques "trou" est superieur ou egal a un, on peut 
d<§finir un critere tel que, lorsqu'en un point de l'objet 
le critere est satisfait, on attribue au point considere 
une resistance mecanique du materiau nulle , de fagon a 
simuler la presence d'une fissure. 
10 L'objet en materiau solide est par exemple une piece 

metallique soumise a une sollicitation mecanique associee 
a une operation de mise en forme par deformation 
plastique, notamment telle que 1 ' emboutissage, l'etirage, 
le laminage, le cintrage ou le forgeage. Cette liste de 
15 procedes de mise en forme par deformation plastique n'est 
pas limitative . 

L'objet peut egalement etre une piece metallique 
soumise & une sollicitation mecanique correspondant a 
l 1 utilisation de la piece. L 1 utilisation de la piece 
20 pouvant engendrer notamment des phenomenes de fluage, de 
fatigue, ou de f atigue-f luage . 

Le proc^de selon I 1 invention peut etre utilise pour 
evaluer l'effet des procedes de mise en forme, tels que 
1 1 emboutissage de toles ou 1 ! extrusion d'un tube, sur les 
25 propri6tes des materiaux. Les resultats obtenus peuvent 
etre utilises par exemple dans des modeles de simulation 
de tests de resistance aux chocs de vehicules 
automobiles . 

Le proc6de selon 1' invention peut egalement etre 
30 utilise pour simuler la fabrication de boites metalliques 
telles que des canettes, obtenues par emboutissage, de 
fac?on a en determiner les propri6tes d'emploi en tenant 
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compte des phenomenes tels que les cornes d' emboutissage 
resultant de 1 ' anisotropie du materiau. 

Le procede selon 1' invention peut aussi etre utilise 
pour simuler la fabrication de tubes d'echangeurs de 
5 chaleur en alliage de nickel pour les generateurs de 
vapeur ou pour simuler la fabrication de gaines en 
alliage de zirconium pour combustible de centrales 
nucleaires, de fagon a mieux determiner les proprietes 
initiales de ces pieces. 
10 L' introduction de briques "joint de grain" et de 

briques "trou" dans le modele permet de simuler les 
endommagements microscopiques du materiau et done de 
prevoir les ruptures ou les risques de rupture en 
service. De ce fait, le procede selon 1' invention peut 
15 egalement etre utilise pour <§tudier le comportement et 
evaluer la duree de vie des installations. En particulier 
dans le domaine des centrales nucleaires pour prevoir la 
duree de vie de composants tels que les cuves, les 
generateurs de vapeur, les tuyauteries, les gaines de 
20 crayons combustibles, tous composants soumis a des 
charges complexes . 

Pour les centrales thermiques classiques, le procede 
peut etre utilise pour prevoir la duree de vie 
d' equipements soumis a des ph^nomdnes de fluage. 
25 Dans 1 1 aeronaut ique, le procede peut etre utilise 

pour la prevision de la duree de vie des disques et des 
ailettes des reacteurs soumis a des phenomenes de fluage 
ou de fatigue-fluage. Le fluage est le mode 
d ' endommagement d'une piece soumise a un chargement 
30 maintenu pendant une longue duree. La fatigue est le mode 
d ' endommagement d'une piece soumise a un grand nombre de 
cycles de chargement successifs. La fatigue-fluage est la 
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combinaison de ces deux modes d' endommagement , par 
exemple lorsque les cycles de chargement comportent un 
maintien de longue duree a contrainte 61evee. 

Ces utilisations ne sont donnees qu'a titre d' exemple 
5 et l'homme du metier saura adapter le procede £ toute 
situation particuli&re faisant intervenir des phenomfenes 
de plasticity ou de fluage. 

L' invention va §tre maintenant decrite plus en detail 
en regard des figures annexees dans lesquelles : 
10 - la figure 1 represente les six syst&mes de 

glissement utilises par le modele polycristallin 
selon 1 1 invention . 

la figure 2 represente les six systemes de 
glissements utilises pour representer le 
15 comportement d'un cristal cubique. 

- la figure 3 est la figure de poles des briques 
"grain" representant un mat6riau isotrope. 

- la figure 4 represente les six systemes de 
glissement utilises pour representer le 

20 comportement d'un cristal hexagonal compact. 

- la figure 5 est la figure de poles des briques 
"joint de grain" d f un module d 1 endommagement 
intergranulaire . 

- la figure 6 represente 1 1 affectation des briques 
25 aux points de Gauss d'un element de calcul par 

elements finis. 

- la figure 7 represente 1 1 affectation des briques 
aux points de Gauss d'un element d'une coque. 

- la figure 8 est la figure de poles des briques 
30 "grain" utilisees pour representer la texture d'un 

mat6riau initialement isotrope. 
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15 



- la figure 9 montre les resultats de calculs 
correspondant a un essai de traction biaxiale sur 
un tube en Zircaloy recristallise . 
Alors que dans les modeles de comportement 
polycristallins connus on utilise autant de systemes de 
glissement qu'il en existe dans le cristal reel, de douze 
a trente selon la nature du cristal, les inventeurs ont 
decouvert de facon nouvelle qu'une bonne representation 
de la deformation d'un cristal pouvait etre obtenue en 
utilisant. seulement six systemes de glissement distincts 
bien choisis. 

Alors que les modeles polycristallins connus 
considerent qu'il est necessaire d'utiliser autant de 
grains qu'il y a d' orientations effectives des grains 
dans le materiau reel, les inventeurs ont decouvert de 
facon nouvelle, qu'il est possible d'obtenir une bonne 
representation du comportement du materiau en utilisant 
un nombre de grains restreint correspondant a au plus dix 
orientations differentes des grains. 

Les inventeurs ont egalement decouvert de facon 
nouvelle qu'il etait possible dans tous les cas de 
chargement, notamment cycliques, et pour les briques 
• "grains" comme pour les briques "joint de grain", 
d'obtenir une bonne representation des materiaux en 
25 utilisant des modeles de changement d'echelle permettant 
de faire le passage de l'echelle macroscopique a 
l'6chelle microscopique et inversement qui ne fasse pas 
intervenir de tenseur supplemental^ mais qui, a la place 
de ces tenseurs, utilise des parametres scalaires. 
30 En outre, les inventeurs ont constate de facon 

nouvelle qu'il n' etait pas indispensable d'af fecter 
1' ensemble des grains du modele de comportement a 
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1' ensemble des points de calcul, et dans le cas du calcul 
par la m§thode des Elements finis, a l 1 ensemble des 
points d ? integration de la methode par elements finis, 
mais qu'il etait possible d'obtenir une precision 
5 excellente en affectant un nombre reduit de grains a 
chacun des points de calculs, et dans le cas de calcul 
par methode des elements finis, un grain a chacun des 
points de calcul de la methode consideree. 

Enfin, les inventeurs ont decouvert de fa<?on 
10 nouvelle, qu'il etait possible de simuler 1 ' endommagement 
& I'echelle microscopique en utilisant des modeles de 
joint de grains simples qui permettent . de simuler les 
phenomenes d' endommagement intergranulaires ou d'utiliser 
des modules de cavites pour modeliser 1 1 endommagement 
15 ductile des materiaux. 

Ainsi, 1' invention se distingue des methodes de 
calcul numerique connues par la nature du modele de 
comportement polycristallin de materiau, et par la fagon 
de l'incorporer dans la methode de calcul numerique. 
20 Af in de bien f aire comprendre 1 1 invention, on va 

considerer le cas particulier de la methode de calcul par 
Elements finis. On va tout d'abord rappeler quelques 
notions sur les techniques de calculs par cette methode 
et sur les modules de comportement de materiaux de type 
25 polycristallin, puis on va expliquer comment, selon 
1 ? invention, on construit un modele de comportement 
polycristallin performant, et on va indiquer comment on 
incorpore ce modele de comportement polycristallin dans 
la methode de calcul par Elements finis. 
30 De fa<?on classique, pour simuler par un calcul par 

Elements finis le comportement d'un objet sounds & des 
solicitations mecaniques, on construit une 
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representation de 1'objet materiel consid§re k 1 1 aide 
d'un maillage, et on identifie les sollicitations 
mecaniques exterieures auxquelles il va §tre soumis. Puis 
de fac?on connue en elle-meme, on calcule par iteration, 

5 les contraintes et les deformations a des instants 
successifs separ^s par un pas de temps At. Ces calculs 
font intervenir un modele de comportement du mater iau. 

Dans le cas d'un calcul par elements finis utilisant 
un modele de comportement du mat^riau du type "modele 

10 polycristallin", a un instant determine t, on connait 
pour chacune des mailles du maillage, le tenseur des 
contraintes macroscopiques S et le tenseur des 

deformations macroscopiques E qu'on peut decomposer en 
un tenseur des deformations Plastiques E e et un tenseur 

15 des deformations plastiques B P et pour chacun des grains 

du modele polycristallin, le tenseur des deformations 

p 

plastiques macroscopiques S g . 

A l ! aide d'un modele de localisation permettant le 
passage de l'6chelle macroscopique a l'echelle 
20 microscopique, on determine le tenseur des contraintes 

dans chacun des grains du modele polycristallin. 

A l'aide d'un modele de comportement anelastique tel 
qu'un modele plastique ou viscoplastique (de tels modeles 
sont connus en eux-memes de l'homme du metier), et des 
25 systemes de glissement du grain, on determine le tenseur 

.p 

des vitesses de deformation microscopiques de chaque 

grain . 
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En integrant les vitesses de deformations sur un 
intervalle de temps At, on calcule les deformations 
plastiques microscopiques £^ & l 1 instant t + At. 

En utilisant un modele d ? homogdneisation permettant 

5 le passage de l'echelle microscopique a l'echelle 

macroscopique, on determine alors les tenseurs des 

deformations E et des contraintes S macroscopiques a 
= — 

l 1 instant t + At. On peut alors recoinmencer le calcul et 
ainsi determiner devolution au cours du temps des 
10 contraintes et des deformations dans l'objet etudie. 

Pour mettre en oeuvre une methode de calcul numerique 
du type de celle qui vient d'etre decrite, on utilise 
trois reperes : 

- un repere global, ou repere du laboratoire, (Li, L 2 , 
15 L 3 ) qui est absolu et qui sert de reference pour la 
description de l'objet a etudier. Selon la nature de 
l'objet <§tudie, ce repere permet de reperer chaque point 
de l'objet par ses coordonnees cartesiennes ou par ses 
coordonnees cylindriques, ou par ses coordonnees 
20 spheriques . 

des reperes locaux macroscopiques orthonormes, 
associes a chacun des points de l'objet etudie. Ainsi, a 
chaque point, et en particulier, & chaque point pour 
lequel on effectue un calcul, on associe un repere local 
25 macroscopique. Ce repere est un repere orthonorme dont la 
position et 1 ' orientation par rapport au repere global 
peut evoluer au fur et £l mesure que l'objet etudie se 
def orme . 

des reperes locaux microscopiques orthonormes 
30 associes a chacune des "briques" du modele polycristallin 
et qui servent a decrire les deformations microscopiques 
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cie chaque brique. L ! orientation du repere local 
microscopique d ! une "brique" associee & un point de 
l 1 ob jet etudie represente 1 1 orientation de ce . grain par 
rapport au repere local macroscopique du point considdre. 
5 Par la suite on appellera (X lf X 2 , X 3 ) un repere local 

macroscopique et (x x , x 2 , x 3 ) un repere local 
microscopique . 

On va maintenant d6crire le modele polycristallin 
utilis6 dans 1' invention. 
10 Selon l f invention, pour construire le modele 

polycristallin, on utilise des briques <§lementaires de 
trois natures differentes : 

des briques "grains" qui representent le 
comportement de grains cristallins ; 
15 - les briques "joints de grains" qui representent le 

comportement de joints de grains ; 

- les briques "trou" qui representent des cavites. 
Les briques "grains" permettent de representer le 
materiau en 1' absence d 1 endommagement microscopique. Les 
20 briques "joints de grains" permettent de representer des 
endommagement s du type intergranulaire . Les briques 
"trou" permettent de representer les endommagement s du 
type ductile. 

On va maintenant deer ire d'abord les briques 
25 "grains", puis les briques "joints de grains" et enfin 
les briques "trous". 

De fagon connue, a chaque brique "grain" qui est 
l 1 equivalent des "grains" utilises dans les modeles 
polycristallins connus, on associe des systemes de 
30 glissements, des lois de comportement du materiau pour 
chacun des systemes de glissement, une orientation de la 
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brique "grain" par rapport a un repere lie a l ! objet que 
1 1 on veut etudier, et enfin une fraction volumique f. 

Selon l 1 invention, et quelle que soit la nature 
cristallographique du materiau solide considere, on 

5 attribue k la brique "grain" six systemes de glissement. 
De fagon connue, chaque syst&me est d6fini par la normale 
a un plan de glissement et par une direction de 
glissement, ce qui permet d'associer ci chaque systeme de 
glissement une matrice d ! orientation. Dans un repere 

10 local microscopique orthonorme (x if x 2 , x 3 ) li£ a la 
brique "grain", selon 1 ! invention, on choisit les 
systemes de glissement de telle sorte que les doubles des 
matrices d 1 orientation, rn^es systemes de glissement ne 

soient constitutes que de 0, 1, -1. 
15 si n s est le vecteur normal au plan du systeme de 

glissement d'indice s, et l s le vecteur de la direction 

de glissement de ce syst&me de glissement, on a par 

definition : 

2^ = ( n s ®L s + L®Hs ) 
20 Avec ces definitions, les systemes de glissement, 

represents a la figure 1, sont les suivants : 
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A la figure 1, on a represents un cube ABCDEFGH dans 
le repere local microscopique (xi, x 2 , x 3 ) . 

Le systeme de glissement s=l, correspond a un 
glissement selon le vecteur li, parallele a l'axe x 2 , dans 



WO 2005/015167 



17 



PCT/FR2004/002012 



un plan tel que le plan CDEH qui est perpendiculaire au 

vecteur ni et parallele au plan (x 2 , x 3 ) . 

Les systemes de glissement s=2 et s=3, sont 

comparables au systeme 3=1, et correspondent a des 
5 glissements paralleles aux axes x 3 et xi, le long de plans 

parall&les aux plans (x 3 , Xi) et (xi, x 2 ) , respectivement. 
Le systeme de glissement s=4 correspond a un 

glissement selon le vecteur 1 4 parallele a la diagonale 

DF du carre ADEF situe dans le plan (xi, x 2 ) , le 
10 glissement se faisant dans un plan tel que le plan DFGC 

perpendiculaire au vecteur n 4 . Ce plan est parallele a 

l'axe x 3 et forme un angle de 45° avec les plans (xi, x 3 ) 

et (x 2 , x 3 ) . 

Les systemes de glissement s=5 et s=6 sont 

15 comparables au systeme s=4 et correspondent a des 

glissements paralleles aux diagonales FB et BD, le long 

des plans FBCE et BDEG, respectivement. Ces plans sont 

paralleles aux axes xi et x 2 respectivement, et forment 

des angles de 45° avec les plans (x 2 , xi) et (x 3/ xi) pour 

20 l f un, et avec les plans (xi, x 2 ) et (x 3 , x 2 ) pour l 1 autre. 

• 

Avec de tels systemes de glissement, si y 

s 

represente la vitesse de glissement dans la direction l s 
du systeme de glissement s, le tenseur des vitesses de 
deformation plastiques de la brique "grain" s'ecrit : 

25 
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Pour determiner les vitesses de glissement y s sur 

chacun des systfemes de glissement s, on utilise un modele 
de comportement qui utilise les relations de 
viscoplasticite classiques decrites par exemple dans La 

5 Revue Europeenne des elements finis Vol 3 - n°4/1994, 
pages 515 a 541. Utilisation de modules polycristallins 
pour le calcul par Elements finis. Georges Cailletaud - 
Philippe Pilvin et qui font intervenir de fac?on connue la 
scission resolue x 3 sur le systeme s. La scission resolue 

10 t s s 1 ecrivant : 

T s =cr : m s 
=z — 

Dans cette equation, cr^ est le tenseur des 

contraintes dans le grain considere. 

Pour modeliser le comportement d'un materiau 

15 determine, il faut alors choisir des briques "grain" du 
modele polycristallin de telle sorte que ces briques 
representent d'une part le comportement d'un grain du 
materiau, compte-tenu de sa nature cristallographique 
reelle et d 1 autre part d'un agregat de tels grains, 

20 compte-tenu de la repartition reelle des orientations des 
grains dans le materiau, c'est-a-dire de la texture du 
materiau. 

Pour faire ce choix, on associe a chaque systeme de 
glissement une scission critique telle que la brique 

25 "grain" obtenue ait un comportement identique a celui du 
cristal reel dont est constitue le materiau, puis on 
choisit le nombre de briques "grains", leurs orientations 
et leurs fractions volumiques de telle sorte que 
1 ! agregat obtenu ait un comportement bien representatif 

30 du materiau reel. 
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Pour repr^senter 1 1 orientation des grains, on 
considere un rep&re local macroscopique (Xi, X 2 , X 3 ) et un 
repere local microscopique (xi, x 2 , x 3 ) lie a chaque 
brique "grain". L 1 orientation d'une brique "grain" est 

5 definie comme etant 1 1 orientation de son repere local 
microscopique (xi, x 2 , x 3 ) par rapport au repere local 
macroscopique (Xi, X 2 , X 3 ) . 

De fa?on connue en elle-meme, on considere que le 
repere local microscopique (xi, x 2 , x 3 ) se deduit du 

10 repere local macroscopique (Xi, X 2 , X 3 ) par trois 
rotations : 

- une rotation d 1 angle \|/ du repere (Xi, X 2 , X 3 ) autour 
de l f axe X 3 qui conduit a un repere (X'i, X f 2 , X 3 ) , 

- une rotation d 1 angle <J> autour de l'axe X f i qui 
15 conduit a un repere (X'i, X" 2 , X f 3 ) 

- une- rotation d 1 angle q> du repdre (X f i, X" 2 , X' 3 ) 
autour de l'axe X ! 3 qui conduit a un repere (X"i, X"' 2 , 
X' 3 ) ; ce dernier repere est le repere local 
microscopique (xi, x 2 , x 3 ) . 

20 Le repere local microscopique (x X/ x 2 , x 3 ) est done 

d6fini par les trois angles de rotation *P, <|>, cp, qui ont 
permis de l'obtenir. Ces angles, sont habituellement 
appeles "angles d'Euler". 

Cette methodologie de choix des briques "grains" va 

25 etre illustree par deux exemples que 1 f homme du metier 
saura transposer & tout cas particulier. 

Le premier exemple, concerne le cas d ! un materiau 
isotrope constitu6 de cristaux de structure cubique 
centree ou faces-centrees tels que les aciers. 

30 Lorsque le materiau est constitu£ de cristaux de 

structure cubique, les systemes de glissements peuvent 
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etre regroup^s en deux sous-families correspondant 
chacune a une scission critique. Ces deux families de 
systdmes de glissement correspondent d'une part aux 
glissements sur les faces du cube, et d 1 autre part aux 

5 glissements sur les plans a 45° de ces faces. 

Ces six syst&mes de glissement adaptes a un cristal 
cubique sont representes a la figure 2. Les systemes s=l, 
s=2 et s=3, correspondent aux glissements sur les faces 
du cube et correspondent a la meme scission critique , 

10 puisque, hormis 1 ' orientation, ces faces sont identiques. 
Les systemes s=4, s=5 et s=6, correspondent aux 
glissements sur les plans k 45° par- rapport aux trois 
directions des aretes du cube. Ces trois systemes 
correspondent a la meme scission critique, puisque, 

15 hormis 1 1 orientation, ces plans sont identiques. 

Afin de bien representer le caractere isotrope du 
materiau, on considere un modele constitue de sept 
briques grains dont les orientations sont choisies de 
fagon £ couvrir au mieux une sphere. 

20 On utilise alors 7 briques "grain" choisies de telle 

sorte que la figure de poles correspondante soit 
symetrique . 

Pour cela, on construit un premier groupe de trois 
briques "grain" et un deuxieme groupe de quatre briques 
25 "grain" en procedant comme suit : 

- pour le premier groupe, on considere une premiere 
brique "grain" primitive d ! indice g=l, obtenue par une 
rotation d 1 angle ¥=71/4 autour de l'axe X 3 . Pour obtenir 
un comportement identique selon les trois axes, on 
30 symetrise la figure de p61es en ajoutant deux briques 
"grain" d'indice g=2 et g=3 obtenues par rotation de 7t/4 
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autour des axes Xi et X 2 respect ivement . Les coordonnees 
des pdles de ces grains sont alors 0, ± 1, ± 1/V2 et 
± (-v/2-1) . 

- pour le deuxieme groupe, on considere une deuxieme 
5 brique "grain" primitive dont les poles sont situes sur 
trois grands cercles passant par un des axes Xi et 
formant un angle de 45° avec les deux autres axes, et 
dont les coordonnees sont ± 0,2, ± 0,4 et ± 0,5. Puis on 
symetrise la figure de poles en ajoutant trois briques 
10 "grains" dont les poles se deduisent des poles de la 
deuxieme brique "grain" primitive par symetrie par 
rapport aux axes Xi, X 2 et X 3 . 

On obtient ainsi sept briques "grains" auxquelles on 
attribue la fraction volumique f = 1/7 et dont les 
15 angles, *F , <)>, et cp, exprimes en radians, sont : 

grains: 12 3 4 5 6 7 

tc/4 0 K /2 0,464 w / 4 2,034 -2, 678 

(j) 0 7t/4 n/ 4 0,841 1,231 0,841 0,841 

cp 0 0 -tc/2 0, 464 -tc/4 -1,107 2,034 



Les sens de deux des trois axes xl, x2 et x3 peuvent 
etre inverses, et les trois axes peuvent etre soumis a 
20 une permutation circulaire sans modifier les resultats du 
modele. Les angles du tableau correspondent a un choix 
particulier des sens des axes et de la permutation. 
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La figure de poles obtenue est representee & la 
figure 3. Sur cette figure, les poles rep^res li, 1 2 / I3 
representent les poles de la premiere brique "grain". Les 
poles rep6r§s 2i, 22, 2 3 representent les poles de la 
5 deuxieme brique "grain" , et ainsi de suite. On voit en 
particulier que les poles 4 3 , 5 2 , 61 et 7 3 ont pour 
coordonnees (+ 0,5 ; + 0,5). Les poles 4i, Az, 5i, 53, 62, 
63, 1\ et 7 2 ont pour coordonnees (+ 0,2 ou + 0,4). 

Ce module est universel pour les materiaux isotropes 
10 a structure cubique centree ou faces-centrees. 

A titre de deuxieme exemple, on considere un materiau 
anisotrope ayant une structure hexagonale compacte. Un 
tel materiau est par exemple du Zircaloy recristallise . 
Pour une telle structure cristallographique hexagonale, 
15 on considere trois groupes de systemes de glissement 
correspondant a trois scissions critiques distinctes. Le 
premier groupe rassemble les systemes 2 et 3 qui 
correspondent au glissement basal, le deuxieme groupe 
rassemble les systemes 1 et 4 qui correspondent au 
20 glissement prismatique, et le troisieme groupe rassemble 
le systeme 5 et 6 qui correspondent aux glissements 
pyramidaux. 

Les syst&mes de glissement destines a representer un 
cristal hexagonal compact sont representes a la figure 4 . 

25 Sur cette figure, les axes xi, x 2 , x 3 constituent le 
repere orthogonal de reference de la brique "grain", 
c ? est-a-dire le repere local microscopique . Le prisme a 
base hexagonale ABCDEFGHI JKL est une structure 
qui dans le cadre du modele a un comportement identique a 

30 celui d'un cristal hexagonal, mais n'a pas la meme 
geom6trie que celle d'un cristal de structure hexagonale. 
Les systemes de glissement : 



WO 2005/015167 



23 



PCT/FR2004/002012 



- s=2 rep6r6 P 2 sur la figure, et s=3, rep6re B 3 
correspondent au glissement basal ; il leur est 
attribue une premiere valeur de scission critique ; 

- s=l repere Pi et s=4 repere P 4 , correspondent au 
5 glissement prismatique ; il leur est attribu6 une 

deuxieme valeur de scission critique ou deux valeurs 
distinctes . 

- s=5 reper<§ 7i 5 et s=6 repere n Sf correspondent aux 
glissements pyramidaux ; il leur est attribue une 

10 valeur de la scission critique particuli^re. 

En outre, il est necessaire de choisir un certain 
nombre de briques "grains" ayant des orientations 
permettant de representer 1 1 anisotropie globale du 
mater iau. Dans le deuxieme exemple, qui concerne en 

15 particulier un tube en Zircaloy recristallise, 
l f invention utilise six briques "grains" ayant chacune 
une fraction volumique de 1/6, dont les orientations par 
rapport au repere local macroscopique (Xi, X 2 , X 3 ) lie au 
point de l'objet considere, sont les suivantes : 

20 
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7t-^i 
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x/2 
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Les orientations de ces six grains dependent de cinq 
parametres qui correspondent a cinq angles d'Euler *¥ lf 
x i*2r <Pi f <P2 et cp 3 . Ces cinq parametres peuvent etre 
25 determines par une methode d 1 identification numerique en 
comparant les resultats de calcul numerique a des 
resultats d'essais de traction-torsion. Ces methodes 
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d ? identification num6rique sont des m6thodes connues de 
l'homme du metier. 

Pour modeliser l 1 evolution de la texture d'un 
materiau polycristallin, c'est-&-dire devolution des 
5 orientations des grains, on consid&re de fagon connue la 
rotation Q du r6seau cristallin d'un grain par rapport a 

un repere dit "corotationnel". La vitesse de rotation du 
reseau cristallin est donn6e par : 

n 

— - 13 = s 

10 ouapparait le tenseur de rotation : 

£ 5 = (l 5 ®n s -n 5 <8>l 5 )l2 

Le repere corotationnel suit la rotation moyenne de 
la matiere R . De fagon connue, sa vitesse de rotation 
par rapport au repere du laboratoire (LI, L2, L3) est 
15 calculee a partir des gradients des vitesses v des points 
materiels du corps solide deformable : 

R:R T = (gradv - grad T v) / 2 

Dans 1 1 exemple du materiau initialement isotrope 
constitue de 7 briques "grain", il peut etre interessant, 

20 pour certaines geometries de pieces et certains 
chargements particuliers coinme le laminage d'une tole, de 
maintenir une symetrie orthotrope apres evolution de la 
texture. Pour cela et de fagon nouvelle, on dissocie les 
grains 1 a 3, obtenus par rotation de 7i/4 autour des 

25 trois axes X3, XI et X2 du rep&re local macroscopique, en 
six grains obtenus par des rotations de n/S (grains 1 a 
3) et 37C/8 (grains 8 a 10) autour de ces axes . De plus, 
les axes x2 et x3 du repere A l'echelle microscopique du 
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grain 8 sont inverses, les axes xl et x3 du grain 9 sont 
inverses et les axes xl et x2 du grain 10 sont inverses, 
ce qui est equivalent a des rotations d f angle n autour 
des axes xl, x2 et x3 respectivement pour les grains 8, 9 
5 et 10. On obtient ainsi dix briques "grains", dont les 
angles d'Euler v|/, (J) et 9 et les fractions volumiques f 
sont : 



grains 


1 
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4 
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10 




7l/8 
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%/2 


0,464 


w/4 


2,034 
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7C/8 
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lbt/8 


37t/8 
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-1,107 
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nil 
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1/14 


1/14 


1/14 


1/7 


1/7 


1/7 


1/7 


1/14 


1/14 


1/14 



10 La figure de poles obtenue dans l'etat initial, avant 
evolution de texture, est representee a la figure 8. Pour 
obtenir une figure de poles orthotrope apres evolution de 
texture, les dix briques "grains" 1 a 10 peuvent etre 
completees par dix briques "grains" obtenues en 

15 substituant les unes aux autres les normales aux plans de 
glissement n s et les directions de glissement [ s des 6 
systemes de glissement, ce qui revient & conserver les 
tenseurs d f orientation m^ f et a inverser les tenseurs de 

rotation q . 

20 On va maintenant decrire les briques "joint de 

grains", destinees a modeliser les endoramagements 
intergranulaires lies k l'ouverture des joints de grains. 
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Pour mod<§liser les endommagements intergranulaires , 
on utilise des briques "joint de grains" dont le 
comportement est represents par au moins un systeme 
d'ouverture et eventuellement au moins un systeme de 
5 glissement. 

Un systeme d'ouverture se distingue d'un systeme de 
glissement par le fait qu f il conduit k un tenseur de 
deformations plastiques dont la trace est positive. 

Comme pour les briques "grains", les systemes 
10 d'ouverture ou de glissement sont definis par des 
matrices d ! orientation m se deduisant des normales aux 

plans de deformation (glissement ou ouverture) « 5 et des 
directions de deformation l s . 

A titre d'exemple, on considere les deux cas suivants 

15 : 

1) la brique "joint de grain" comport e un systeme 
d'ouverture (s=l) et deux systemes de glissement 
intergranulaire (s=2 et s=3) . Ces systemes sont definis 
de la fagon suivante : 

20 

s = 1 2 3 

n s « (1 0 0) (1 0 0) (1 0 0) 



Is = 


(1 


0 


0) 




f 2 


0 


0 


2 m a » 


0 


0 


0 






0 





(0 

(5 
i 
o 

vl 



l 

0 
0 



0) 

0 ^ 
0 
0 



(0 



ro 

0 

1 

V. 



0 
0 
0 



1) 

1^ 

0 
0 



WO 2005/015167 



27 



PCT/FR2004/002012 



2) La brique "joint de grain" comprend trois systemes 
d'ouverture (s=l, s=2, s=3) et trois systemes de 
glissement intergranulaire (s=4, s=5, s=6) : 
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Le choix des briques "joint cle grain" depend de la 
nature des sollicitations auxquelles est soumis l'objet 
dont on veut etudier le comportement . 

10 Lorsque l'objet dont on veut etudier le comportement 

est soumis & des sollicitations dans une seule direction, 
XI, on peut considerer une seule brique "joint de grain" 
perpendiculaire a la direction XI, dont le syst£me 
d'ouverture correspond £ la direction XI et les deux 

15 autres systemes de glissement correspondent a X2 et X3 . 
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Pour un chargement dans une direction quelconque, on 
peut soit utiliser trois briques "joint de grain" 
perpendiculaires a chacune des directions XI, X2, X3, 
soit utiliser une seule brique "joint de grain" 

5 comportant trois systemes d'ouverture et trois systemes 
de glissement. 

Enfin pour la description fine de 1 1 endommagement par 
une sollicitation multiaxiale, on peut considerer sept 
briques "joint de grain" ayant chacune un syst^me 

10 d'ouverture et deux systemes de glissement, et une 
orientation definie par les deux angles d'Euler de la 
normale n du systeme d'ouverture, 

Les orientations des directions d'ouverture des 
briques "joint de grain" peuvent etre representees par 

15 une figure de poles reproduite a la figure 5, dans 
laquelle les sept orientations sont reper6es 1, 2, 3, 4, 
5, 6 et 7. 

Les fractions volumiques des briques "joint de grain" 
sont choisies de telle sorte que leur fraction volumique 

20 totale soit egale a six fois le rapport de 1'epaisseur 
mecanique du joint de grain & son diametre. A titre 
indicatif, cela represente quelques pour-cent. 

Pour decrire le comport ement des joints de grains, on 
utilise des modeles visco-plastiques classiques qui 

25 permettent de determiner les vitesses d'ouverture des 
joints de grain. Ces modeles sont par exemple les modeles 
decrits par Onck et Van der Giessen dans Journal of the 
Mechanics and Physics of Solids-Vol 47, 1999, pages 99- 
139-"Growth of an initially sharp crack by grain boundary 

30 cavitation" . 

On va maintenant decrire une brique "trou", destinee 
a modeliser 1 ' endommagement ductile. Une telle brique 



WO 2005/015167 



29 



PCT/FR2004/002012 



"trou" est caracterisee par une vitesse de deformation 
plastique purement volumique representee par les 
equations : 



el =l/3D ie xp 



m 



E p 



p = exp(-32?£) 



Em "■ f 



f = l-(l-f 0 )exp (-3E*,) 



dans lesquelles : 

Di et Gi sont des coefficients scalaires caracteristiques 
du materiau , 

10 s p m est la vitesse de deformation plastique volumique 
microscopique . 

Sm est la contrainte macroscopique moyenne, <§gale au 1/3 
de la trace du tenseur des contraintes. 

Ef g est la vitesse de deformation plastique macroscopique 
15 equivalente, au sens de von Mises. 

E p est la vitesse de deformation plastique macroscopique 

171 

moyenne . 

f 0 est la fraction volumique initiale de la brique "trou" 
f est la fraction volumique de la brique "trou" a 
20 l f instant t. 

p est la masse volumique relative. 



L 
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Cette brique "trou" a une fraction volumique f qui 
est une variable. Quand on ajoute une brique "trou" a un 
modele constitu§ de briques "grain" de fraction volumique 
f g chacune, on consider e que la fraction volumique de 
5 chaque brique "grain" devient (1-f) f g . 

Les differentes briques "grain", "joint de grain" et 
"trou" qui viennent d'etre decrites, permettent de 
modeliser le comportement microscopique du materiau. 
Mais, comme on l f a vu precedemment pour modeliser le 
10 comportement d'un objet, il est necessaire de pouvoir 
passer du niveau microscopique au niveau macroscopique et 
reciproquement . Pour cela, on utilise des modeles appeles 
"modules de changement d ! §chelle" dont le principe 
general est connu de l'homme du metier, mais qui, dans le 
15 cadre de la present e invention, comportent des 
particular ites qu'on va decrire maintenant. 

On considere trois situations selon que le modele ne 
comporte que des briques "grain", que le modele comporte 
des briques "grain" et des briques "joints de grain" ou 
20 que le modele comporte des briques "grains" et une brique 
"trou". 

Lorsque le modele ne comporte que des briques 
"grains", les relations permettant le passage de 
l'6chelle macroscopique & l f 6chelle microscopique et 
25 inversement sont respectivement : 

a) pour le passage de l'echelle macroscopique a 
I'echelle microscopique : 

a g = E + 2 n (l-p)a(E P -/ ) 
avec : 
30 1 = 1 + D(E P )7 
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et 

2 4-5v 



Dans ces equations, 

5 a est le tenseur des contraintes microscopiques dans 

=g 

la brique "grain" d'indice g. 

Z est le tenseur des contraintes macroscopiques • 

est le tenseur des deformations plastiques 

macroscopiques . 

10 / est le tenseur des deformations microscopiques 

dans la brique "grain" d'indice g. 

|x est le module de cisaillement du materiau. 
v est le coefficient de Poisson du materiau. 
D est un parametre propre au materiau, ajustable a 
15 partir d'essais. 

/ p\-" egt la valeur maximale du deuxieme invariant de 

von Mises du tenseur des deformations plastiques 
macroscopiques, atteinte anterieurement . 

b) Pour le passage de l'echelle microscopique a 
20 l f echelle macroscopique : 

N 



N p 

e = j(fi|; 



) 
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dans ces equations , 

f g est la fraction volumique de la brique "grain" 
d'indice g, 

E fl est le tenseur des deformations 61astiques 
macroscopiques . 



10 



15 



Lorsque le modele comporte des briques "grains" et 
des briques "joint de grain", les relations de passage du 
niveau macroscopique au niveau microscopique et du 
passage du niveau microscopique au niveau macroscopique 
sont modifiees de la fagon suivante : 



avec 

i=i+D<E p r 

(X. — eq 

et 



1=2+5D (E p ) 



max 

> 

eq 



Enfin, lorsque le modele comporte des briques 
"grains" et une brique "trou", les relations de passage 
20 du niveau macroscopique au niveau microscopique et 
inversement sont les suivantes : 

a) pour le passage de l'echelle macroscopique a 
I'echelle microscopique : 



25 



et 



S + D x f cr, exp 



pa, 



-eq 
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a Q est le tenseur des contraintes microscopiques dans 

la brique "trou" . 

Di et Ci sont des parametres propres au materiau. 
b) et pour le passage de l'6chelle microscopique a 
5 l'6chelle macroscopique : 



Les briques telles qu'elles ont ete definies plus 
haut, associees aux modeles de localisation qui viennent 
d'etre deer its, constituent le modele polycristallin. Ce 
modele permet de calculer l'effet de sollicitations 

15 mecaniques sur un objet, mais pour cela il doit etre 
integre a un module de calcul par elements finis. 

On va maintenant decrire 1 1 integration du modele 
polycristallin de comportement du materiau dans un modele 
de calcul par Elements finis. 

20 On va tout d'abord decrire le cas d ! un calcul sans 

endommagement , e'est-a-dire dans lequel le modele de 
comportement ne comporte que des briques "grains". Dans 
un tel modele, le nombre de briques "grains" etant 
faible, on affecte une seule brique "grain" a chaque 

25 point de discretisation du calcul de structure par 
elements finis. Les briques "grains" sont affectees aux 
points de Gauss, qui sont les points d 1 integration dans 
le calcul par elements finis. En chacun de ces points, le 



n 



S = f<r n + (1-f) 2 (fg<7 ) 




8" = 2(rf)Yl 

=8 s r 
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nombre de briques "grains" etant egal a un, seules les 
relations scalaires du comportement des briques "grains" 
sont utilis6es. De plus, les vitesses de deformation 
plastique macroscopiques sont confondues avec celles des 
5 briques "grain", II en est de meme de la contrainte 
macroscopique qui est 6gale k la contrainte dans la 
brique "grain" et de la deformation macroscopique qui est 
egale a la deformation dans la brique "grain". II en 
resulte que les ■ relations de localisation et 
10 d f homogeneisation du modele de changement d'echelle 
deviennent inutiles . 

Afin de bien representer la structure du materiau et 
en particulier son homogeneite dans I'espace, les briques 
"grain" sont affectees aux points d 1 integration du 
15 maillage de la structure au moyen d'un algorithme 
garantissant une permutation au sein d'un meme element 
fini et entre les elements finis voisins. Ainsi, les 
briques "grain" sont reparties de fa<?on analogue aux 
grains du materiau reel. Un exemple de repartition est 
20 donne a la figure 6 correspondant au cas d'un materiau a 
structure cubique tel que decrit precedemment . Cette 
figure represente un element fini repere generalement par 
10 comportant six faces en forme de quadrilatere . Les 
sommets et les milieux des aretes correspondent a des 
25 noeuds du modele reperes 11, ce qui fait vingt nceuds. Les 
points d 1 integration ou points de Gauss, reperes 12, sont 
situes a I'interieur de l 1 element fini. Les briques 
"grain" correspondant aux sept orientations sont 
reparties comme I'indiquent les chiffres 1 a 7 associes 
30 aux points de Gauss. Dans le cas ou l f objet etudie est 
une plaque ou une coque on represente cet objet par 
plusieurs couches d'epaisseur identique avec trois points 
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d 1 integration en milieu et en limite de chaque couche, ce 
qui conduit a un nombre de points d 1 integration, egal a 2 
p + 1. Comme on utilise de trois a cinq couches le nombre 
de points d' integration est de sept, neuf ou onze. Ainsi, 
5 la plaque ou la coque comporte toutes les briques "grain" 
du modele dans son epaisseur, comme un materiau reel. Une 
permutation garantit une bonne representation du. materiau 
reel dans les directions tangentes k la plaque ou a la 
coque. Une permutation garantit egalement une bonne 
10 representation moyenne du materiau reel dans 1 ! epaisseur 
de la plaque ou de la coque. La repartition des briques 
"grain" est representee a la figure 7. L ! Element fini 20 , 
comporte trois couches 21, 22 et 23. Les points 
d f integration sont regroupes en quatre ensembles 24, 25 , 
15 26 et 27 comportant chacun sept briques "grain" dont les 
orientations sont reperees de 1 a 7 . Dans certains cas il 
peut etre interessant d'affecter plusieurs briques 
"grain" a un des points d 1 integration voire toutes les 
briques "grain" du modele a chaque point d ! integration. 
20 Dans ce cas, les calculs sont un peu plus longs. 

En particulier, dans le calcul par elements finis 
avec evolution de texture utilisant les dix briques 
"grain" precedemment definies pour un materiau 
initialement isotrope, on affecte a certains points 
25 d 1 integration deux briques "grains" correspondant aux 
grains dissocies 1 et 8, 2 et 9, 3 et 10. Si l'on utilise 
les dix briques "grains" supplementaires obtenues en 
substituant les normales aux plans de glissement n s et 
les directions de glissement lg , on affecte ces briques 
30 aux memes points d 1 integration que les briques 1 k 10. Le 
nombre de comportements distincts des points 
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d' integration reste ainsi limite k sept, comme dans le 
cas d'un calcul sans evolution de texture pour un 
mat^riau isotrope. 

Lorsqu'on considere le cas dans lequel les joints 

5 peuvent etre endommages, c'est-a-dire lorsqu'on utilise 
un modele comportant des briques "joint de grain", les 
fractions volumiques des briques "joint de grain" etant 
tres faibles, pour obtenir une bonne representation du 
materiau on associe une brique "grain" et une brique 

10 "joint de grain" k chaque point d' integration du modele 
de calcul par elements finis, et on utilise un algorithme 
de permutation de telle sorte que les couples brique 
"grain" /brique "joint de grain" ne soient pas tou jours 
les memes. Dans ce cas, les relations de changement 

15 d'echelle restent necessaires mais sont tres simplifiees. 
Ces relations sont : 



20 



fg + f j = 1 

a = Z + 2\x (1-P) [ti(1 ® 1) + a l]fj(gj-£j) 
a.= | + 2[i (1-P)[ti(1 ® 1) + a J]f g (5j-5j) 



25 Dans le cas ou I 1 on veut modeliser 1 ' endommagement 

ductile, c'est-a-dire lorsqu'on utilise un modele 
comportant des briques "grain" et une brique "trou", on 
associe une brique "grain" et une brique "trou" k chacun 
des points d' integration . Dans ce cas egalement, les 
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relations de localisation sont tres simplifies et ces 
relations sont : 

Erfev = (1-f ) £ g 



15 



20 



25 



1- = f £ 0m 



<r = — — 2 + DipCTi exp! 



fill 

V. per, J S„ 



Avec ce modele, il est alors possible de calculer 
Involution entre un instant 0 et un instant t, des 
proprietes du mater iau dont est const itue l'objet pour 
10 lequel on fait 1' etude des contraintes et des 
deformations a l'etat macroscopique . Cependant, pour 
faire un tel calcul, il est necessaire de pouvoir 
determiner les differents parametres qui interviennent 

dans le modele . 

Pour determiner ces parametres, une methode 
d' identification connue en elle-meme, selon laquelle on 
commence d'abord par rassembler dans une base de donnee 
relative au materiau des resultats d'essais de traction, 
d'essais traction/torsion ou encore d'autres essais 
mecaniques tela que des essais de fluage. Puis on simule 
par le calcul les differents essais de la base de donnee 
en utilisant des valeurs a priori des parametres. Puis on 
compare les resultats de ces simulations aux essais 
proprement dits et, en fonction de l'ecart entre les 
resultats experimentaux et les resultats calcules, on 
deduit de nouvelles valeurs pour les parametres. Enfin, 
on reproduit ce processus jusqu'a ce que les ecarts entre 



WO 2005/015167 



38 



PCT/FR2004/002012 



les resultats experimentaux et les resultats simules 
soient infSrieurs a une valeur definie a l'avance. 

A tit re d 1 exemple, on a applique le modele selon 
I 1 invention au cas d f un tube en Zircaloy recristallis6 
5 soumis a un effort de traction biaxee. Dans le repere du 
Laboratoire (Li, L 2 , L 3 ) , on a utilise des coordonn^es 
cylindriques . 

Le modele a 6t6 identifie sur la base de deux essais 
de traction dans les orientations axiales zz et 
10 tangentielle 99 du tube et d'un essai de torsion z9 du 
tube. Le resultat de 1 ' identification est represents a la 
figure 9. Dans cette figure, les triangles reperes 31 
correspondent a des points experimentaux, la courbe 32 
represente les resultats tels qu'ils sont obtenus en 
15 appliquant les criteres de von Mises avec un module 
macroscopique classique de comportement du materiau, la 
courbe 33 represente les resultats obtenus par un modele 
polycristallin selon l f art anterieur, et la courbe 34 
represente les resultats obtenus avec le modele selon 
20 l f invention, Ces courbes montrent que le procede selon 
l 1 invention conduit a une tres bonne concordance des 
resultats experimentaux et des resultats calcules pour 
des tractions dans toutes les directions du plan zz-99, 
ce qui n'est pas le cas de 1 1 application simple des 
25 criteres de von Mises, ni le cas de 1 1 utilisation d'un 
modele polycristallin selon 1 1 art anterieur. 
L ! identification a demande 2 heures de calcul sur une 
station de calcul standard. L 1 identification du modele 
polycristallin selon l ! art anterieur est decrite dans 
30 "Philippe Geyer, Comportement elasto-viscoplastique de 
tubes en Zircaloy-4 : approche expSrimentale et 
modelisation micromecanique, these de l'Ecole des Mines 
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de Paris, 9 decembre 1999". Cette identification a 
necessite plusieurs jours de calcul sur un ensemble de 5 
a 10 stations de calcul travaillant simultanement . Le 
resultat est moins bon qu'avec le modele selon 
5 1' invention non pas parce que le modele polycristallin 
selon l'art anterieur est moins bon que le modele selon 
1' invention, mais parce qu'il conduit a des calculs tres 
longs qui n'ont pas permis d'optimiser 1 ' identification 
des parametres. Cet exemple montre l'interet du procede 
10 selon 1' invention qui reduit de facon considerable les 
temps de calcul et ainsi permet de determiner l'effet de 
solicitations mecaniques sur des objets dans des 
conditions dans lesquelles les methodes selon l'art 
anterieur ne sont pas applicables. 
15 Le procede qui vient d'etre deer it permet de calculer 

1' evolution des proprietes du materiau ainsi que ses 
endommagements . II fait appel a des modeles visco- 
plastiques, mais l'homme du metier pourra utiliser tout 
modele de comportement anelastique compatible avec un 
20 modele polycristallin conforme a 1' invention et 
notamment, les comportement s en fluage et en fatigue 
f luage . 

Les endommagements correspondent soit a 1 ' ouverture 
des briques "joint de grain" soit au grossissement des 

25 briques "trou" selon qu'on s'interesse a un endommagement 
intergranulaire ou a un endommagement ductile. Dans les 
deux cas, l'homme de metier peut definir un critere de 
rupture defini, par exemple, par un seuil d ' endommagement 
tel, qu'au dela, il y a rupture. On peut ainsi prevoir 

30 les conditions de rupture en service. Dans le cas de 
1 ' endommagement intergranulaire, dans le cadre d'une 
methode par elements finis, on peut definir un critere 
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d ! ouverture d'une brique "joint de grain" tel que, 
lorsque ce critere est satisfait, il y a fissure et, au 
point correspondant , la resistance mecanique est nulle. 
Dans le cas de 1 ' endommagement ductile, dans le cadre 

5 d'une m6thode par elements finis, on peut definir une 
fraction volumique critique d'une brique "trou" telle 
que, lorsque ce critere est satisfait, il y a fissure et, 
au point correspondant, la resistance mecanique est 
nulle. Dans les calculs, on consid&re le point auquel ce 

10 critere est atteint en premier. Puis en prenant en compte 
les points voisins sur lesquels les contraintes vont se 
reporter, on peut suivre la croissance de la fissure et 
ainsi determiner les conditions d' apparition d'une 
rupture . 

15 Le proced6 selon 1 ' invention peut etre utilise pour 

determiner 1' incidence du procede de fabrication d'un 
objet comportant une operation de mise en forme par 
deformation plastique a froid, tels que 1 ' emboutissage, 
et en particulier pour determiner les proprietes 

20 mecaniques de 1' objet en chacun de ses points. En effet, 
1' incidence d'une operation de mise en forme par 
deformation plastique sur les proprietes mfecaniques 
depend de la deformation ef f ectivement realisee et cette 
deformation depend en general du point considere. Les 

25 proprietes mecaniques aprhs mise en forme ne sont done 
pas, en general, identiques en tous points de 1' objet. On 
peut alors utiliser 1' evaluation de la repartition des 
proprietes mecaniques de l 1 objet obtenue par le procede 
selon 1' invention, pour determiner le comportement en 

30 service de cet objet. L'avantage du procede selon 
1' invention est que cette evaluation de la repartition 
des proprietes mecaniques d'un objet peut etre faite dans 



WO 2005/015167 



41 



PCT/FR2004/002012 



des conditions tres diverses a partir d'un nombre 
restreint d'essais mecaniques. En effet, pour evaluer les 
proprieties mecaniques d'objets constitues du meme 
materiau, mais realises par mise en forme selon des 
5 modalites diverses, il suffit de realiser les essais 
mecaniques necessaires a la construction du modele 
polycristallin selon 1' invention, representatif du 
materiau, puis appliquer le procede selon 1' invention a 
chacun des cas particuliers envisages. 
10 A titre d'exemple, on peut utiliser le procede selon 

1' invention pour evaluer le comportement au choc d'une 
piece de vehicules automobiles obtenue par emboutissage 
ou par pliage ou par tout autre procede de mise en forme 
induisant des deformations plastiques . Pour cela, on 
15 construit le modele polycristallin representatif du 
materiau dont est constitute la piece. On cale le modele 
polycristallin & l'aide des essais mecaniques de traction 
ou traction/torsion ou cisaillement necessaires. Puis, a 
l'aide d ? une methode de calcul telle qu'une methode de 
20 calcul par elements finis, incorporant le modele 
polycristallin representatif du materiau, on simule 
l f operation de mise en forme et on determine les 
proprietes mecaniques de la piece en chacun de ses 
points. A partir de ce resultat, et en utilisant une 
25 methode de calcul adaptee connue de 1 1 homme du metier, on 
determine, par exemple, les effets sur les pieces d'un 
choc resultant d'une collision. 

Lorsqu'on souhaite determiner le comportement en 
service de differentes pieces realisees dans le meme 
30 materiau, on utilise le meme modele polycristallin 
representatif du materiau pour ef fectuer les calculs de 
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simulation de comportement . Ainsi, il n f est pas utile de 
faire des essais mecaniques pour chaque pi£ce. 

Le procedt selon l 1 invention peut egalement etre 
utilise pour determiner le comportement d'une piece 
5 constitute d'un materiau pouvant presenter une 
transformation de phrase, c'est-a-dire un changement de 
structure cristallographique sous l ! effet des 
sollicitations therirtomecaniques qu'il subit dans les 
conditions de service ou de mise en forme considerees. 
10 C'est le cas, par exemple, d'un alliage a memoire de 

forme isotrope qui possede deux phases, I 1 une stable a 
haute temperature (1 1 austenite) , l 1 autre metastable a 
basse temperature (la martensite) . La diminution de la 
temperature et/ou 1 1 application d'une sollicitation 
15 mecanique induit un passage progressif de l f austenite 
vers la martensite. La transformation inverse est 
possible de fa<?on reversible. Cette transformation de 
phases s ' accompagne d'une deformation anelastique qui 
s'ajoute a la deformation propre de chaque phase prise 
20 independamment . 

La simulation de cette transformation de phases est 
realisee en utilisant des briques "grains" (7 en 
1' occurrence pour un materiau isotrope), chacune etant 
orientee dans le repere local macroscopique (XI, X2, X3) 
25 par ses angles d'Euler § et cp. Le comportement de 

chaque brique est simule par les six systemes de 
glissement precedemment decrits et leurs lois 
d' ecrouissage connues par ailleurs. Les parametres de ces 
lois de comportement correspondent dans un premier temps 
30 a ceux de 1' austenite. 
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En cours de chargement thermomecanique, en fonction 
d ! un crit^re bien defini et connu, par exemple par 
"Leclercq S., Lexcellent C, A general macroscopic 
description of thermomechanical behaviour of shape memory 
5 alloys, J. Mech. Phys. Solids, Volume 44, pp 953-980, 
1996", l'austenite est susceptible de se transformer en 
martensite, k une vitesse donn6e . A chaque etape du 
calcul, on v6rifie si ce crit£re est satisfait, et on 
attribue aux briques "grain" pour lesquelles le critere 
10 est satisfait, les caracteristiques mecaniques de la 
martensite (evolution des parametres des lois de 
comportement des systemes de glissement de la martensite, 
connue par ailleurs) . La deformation de transformation de 
phase qui resulte de cette Evolution, apparait alors 
15 lorsqu'on applique les regies de transition d'echelle du 
modele. 

Le module qui vient d'etre decrit est essentiellement 
un module numerique qui peut etre mis en ceuvre par tout 
moyen approprie, en particulier, par ordinateur en 

20 I 1 integrant dans un code de calcul par elements finis 
connu en lui-meme . 

Ce modele peut etre utilise de la meme fagon dans 
toute methode de calcul de deformation de structure et, 
notamment, dans les methodes de calculs aux differences 

25 finis ou dans des methodes qui ne font intervenir q ! un 
point de I'objet etudie comme c'est le cas dans la 
modelisation du laminage de produits plats. L'homme du 
metier sait faire les adaptations necessaires qui ne 
portent que sur 1 1 integration du modele dans la methode 

30 de calcul. 

On a decrit les exemples d 1 application a des 
materiaux ayant des structures cristallographiques 
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cubiques ou hexagonales. Mais l f invention ne se limite 
pas a ces exemples et l'homme du metier saura I'appliquer 
a d'autres structures cristallines et a des textures 
diver ses. Elle s 1 applique egalement & des. structures 

5 mixtes constitutes de cristaux de natures differentes. 
Pour cela, il suffit de prendre en compte des briques 
"grain" dont une partie correspond a une structure 
cristalline et une autre partie correspond a une autre 
structure cristalline . 

10 En resume, le proctde selon l f invention, applique a 

un objet constitue d'un materiau determine et soumis 
pendant une periode de temps determinee a des 
sollicitations mecaniques et eventuellement thermiques 
determinees, consiste a : 

15 - choisir un modele de comportement du type 

polycristallin conforme a 1' invention, adapte au 
materiau, 

- ajuster les parametres du modele de comportement a 
partir d'un nombre restreint d'essais mecaniques 

20 convenablement choisis, 

- identifier l'effet des sollicitations mecaniques et 
eventuellement thermiques auxquelles l 1 objet est soumis 
pendant la periode de temps consideree, 

- et calculer l'effet des sollicitations en chaque 
25 point de l 1 objet en utilisant une methode de calcul 

incorporant le modele de comportement du materiau. 

Le procede peut etre applique a tout type de materiau 
dont le comportement peut etre represents par un modele 
du type polycristallin, et en particulier a tout materiau 
30 metallique. 

Les sollicitations auxquelles sont soumises l 1 objet, 
qui peuvent etre des efforts ou des deplacements repartis 
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dans l ! espace et dans le temps et qui peuvent etre 
completees par des sollicitations thermiques, peuvent 
etre determinees par tout moyen et en particulier par des 
enregistrements de mesures. 
5 Les m6thodes de calcul utilisees par le procdde sont des 
methodes de calcul numerique mises en ceuvre par 
ordinateur a l'aide de logiciels adaptes. 
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RE VEN D I CAT I ON S 

1. Proc6de pour determiner les effets d'une 
sollicitation m6canique sur un objet constitue d' un 

5 materiau solide selon lequel on calcule les deformations 
et les contraintes engendrees par la sollicitation 
mecanique en une pluralite de points de 1' objet au moyen 
d' une m£thode de calcul numerique en mode anelastique 
dans laquelle le comportement du materiau solide est 

10 represents par un modele de comportement microscopique 
polycristallin utilisant une pluralite de briques 
« grain » dont les deformations sont determinees a partir 
d'une plurality de systemes de glissement propres au 
materiau solide, la trace du tenseur des deformations 

15 microscopiques etant nulle, chaque brique "grain" ayant 
une orientation et une fraction volumique, caracterise en 
ce que : 

- le nombre de briques « grains » est inferieur ou 
egal a dix, 

20 - le nombre de systemes de glissement est inferieur 

ou egal a six, 

en outre, le modele microscopique utilise eventuellement : 

jusqu'a sept briques « joint de grain » ayant 
chacune une orientation et une fraction volumique, et 
25 dont les deformations sont d6terminees a partir d' au 
moins un systeme d' ouverture, la trace du tenseur des 
deformations plastiques microscopiques de chaque brique 
« joint de grain » etant positive, 

- au plus une brique « trou » ayant une fraction 
30 volumique variable, dont les deformations plastiques sont 

purement volumiques et dont la trace du tenseur des 
deformations plastiques microscopiques est positive. 
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2. Procede selon la revendication 1, caracterise en 
ce que les six systemes de glissement associes a une 
brique « grain », definis dans un repere orthonorme 

5 associe a la brique « grain » par les normales aux plans 
de glissement et par les directions de glissement, sont 
tels que le double de la matrice d' orientation associee a 
chaque syst^me de glissement ne soit constitue que de 0, 
1 ou -1. 

10 

3. Procede selon la revendication 1 ou la 
revendication 2, caracterise en ce que pour un materiau 
isotrope k structure cubique, les deformations sont 
determinees a partir de deux families de trois systemes 

15 de glissement, une premiere famille correspondant aux 
faces du cube, et la deuxieme famille correspondant aux 
plans a 45° de ces faces, et en ce que le nombre de 
briques « grains » est egal a sept, les briques 
« grains » etant obtenues par symetrisation de deux 

20 briques « grains » primitives situees sur les grands 
cercles a 45° de la figure de poles. 

4. Procede selon l'une quelconque des revendications 
1 ci 3 caracterise en ce que on determine en outre 

25 l 1 evolution de la texture du materiau en calculant la 
rotation Q du reseau cristallin de chaque brique "grain" 

par rapport a un repere corotationnel . 

5. Procede selon la revendication 4 caracterise en ce 
30 que, pour un materiau isotrope a structure cubique, le 

nombre de briques "grain" est egal a dix. 
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6. Proc6d6 selon I'une quelconque des revendication 1 
a 3, caracterise en ce que on determine en outre 
1 • avancement d'une transformation de phase. 

5 

7. Procede selon la revendication 6, caracterise en 
ce que on attribue a chaque brique "grain", des 
caracteristiques fonction de 1 ' avancement de la 
transformation de phase, 

10 

8. Proced6 selon l'une quelconque des revendications 
1 a 7, caracterise en ce que les orientations et les 
fractions volumiques des briques « grains » sont 
determinees par ajustement a partir d'essais de traction 

15 biaxiale et de traction torsion. 

9. Procede selon l'une quelconque des revendications 
1 a 8, caracterise en ce que, en 1' absence de brique 
« trou », le modele de passage du niveau macroscopique au 

20 niveau microscopique pour une brique « grain », 
permettant de calculer le tenseur des contraintes 
microscopiques a g dans la brique « grain » en fonction du 

tenseur des contraintes macroscopiques s , du tenseur § P 
des deformations plastiques microscopiques dans la brique 
25 « grain » et du tenseur E P des deformations plastiques 
macroscopiques, s'ecrit : 

£ g = 2 + 2yi(l - (5) [ n (l ® 1 ) + a J] (E P - 8 P ) 

avec, l/a = l + Z)^ ax , et Utj = 2 + 5E™, E%* est la valeur 

maximale atteinte par le deuxieme invariant de von Mises 
30 lors de la deformation. 
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10. Proced6 selon l'une quelconque des revendications 
1 a 9, caracterise en ce que la methode de calcul 
numerique est une methode de calcul par elements finis et 

5 en ce que a chaque point d' integration de la methode de 
calcul par elements finis, on associe une seule brique 
« grain » et eventuellement une brique « joint de grain » 
ou une brique « trou », de telle sorte que la repartition 
des briques soit homogene. 

10 

11. Procede selon l'une quelconque des revendications 
1 a 10, caracterise en ce que la methode de calcul 
numerique est une methode de calcul mecanique ou 
thermomecanique, f aisant intervenir le temps . 

15 

12. Procede selon la revendication 11, caracterise en 
ce que on utilise une brique « trou » et qu'a 1' instant 
initial du calcul numerique, la fraction volumique de la 
brique « trou » n'est pas nulle. 

20 

13. Procede selon l'une quelconque des revendications 
1 a 11, caracterise en ce que le nombre de briques "joint 
de grain" est sup^rieur ou egal a un, et en ce que on 
definit un critere tel que, lorsqu'en un point de 

25 l'objet, le critere est satisfait, on attribue au point 
considere une resistance mecanique du materiau nulle, de 
fa<?on a simuler la presence d'une fissure. 

14. Procede selon l'une quelconque des revendications 
30 1 a 8 et 10 a 12, caracterise en ce que on utilise une 

brique "trou", et en ce que on d6finit un critere tel 
que, lorsqu'en un point de 1'objet, le critere est 
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satisfait, on attribue au point considere une resistance 
mecanique du mater iau nulle, de fagon a simuler la 
presence d'une fissure. 

5 15, Procede selon l'une quelconque des revendications 

1 a 14, caracteris6 en ce que l'objet en matdriau solide 
est une piece metallique et en ce que la sollicitation 
mecanique est la sollicitation mecanique associee a une 
operation de mise en forme par deformation plastique 

10 telle que notamment 1' emboutissage, l'etirage, le 
laminage, le cintrage ou le forgeage. 

16. Procede selon 1'une quelconque des revendications 
1 a 15, caracterise en ce que l'objet en materiau solide 

15 est urie piece metallique et en ce que la sollicitation 
mecanique est une sollicitation mecanique correspondant a 
1' utilisation de la pi£ce metallique. 

17. Procede selon la revendication 16, caracterise en 
20 ce que 1' utilisation en service de la piece engendre au 

moins un phenom^ne, pris parmi les phenomenes de fluage, 
de fatigue et de f atigue-f luage . 
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